Mapping the Asymmetric Thick Disk I. A Search for Triaxiality 

Jeffrey A. Larsen 

Physics Department, United States Naval Academy, Annapolis, MD 21402 

larsenOusna . edu 
Juan E. Cabanela 

Department of Physics and Astronomy, Minnesota State University Moorhead, Moorhead 

MN, 56563 

cabanela@mnstate . edu 

Roberta M. Humphreys 
Astronomy Department, University of Minnesota, Minneapolis MN, 55455 

roberta@umn . edu 

and 

Aaron P. Haviland 

Physics Department, United States Naval Academy, Annapolis, MD 21402 



Received 



accepted 



- 2 - 



ABSTRACT 



A significant asymmetry in the distribution of faint blue stars in the inner 
Galaxy, Quadrant 1 (/ = 20° — 45°) compared to Quadrant 4 was first reported 
by Larsen & Humphreys (1996). Parker et al (2003, 2004) greatly expanded 
the survey to determine its spatial extent and shape and the kinematics of the 
affected stars. This excess in the star counts was subsequently confirmed by Juric 
et al. (2008) using SDSS data. Possible explanations for the asymmetry include 
a merger remnant, a triaxial Thick Disk, and a possible interaction with the bar 
in the Disk. In this paper we describe our program of wide field photometry 
to map the asymmetry to fainter magnitudes and therefore larger distances. To 
search for the signature of triaxiality, we extended our survey to higher Galactic 
longitudes. We find no evidence for an excess of faint blue stars at / > 55° 
including the faintest magnitude interval. The asymmetry and star count excess 
in Quadrant 1 is thus not due to a triaxial Thick Disk. 

Subject headings: Galaxy: structure, Galaxy: kinematics and dynamics 
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1. Introduction: The Asymmetric Thick Disk 

Studies of both stars and gas are revealing significant structure and asymmetries in 
their motions and spatial distributions including the bar of stars and gas in the Galactic 
bulge (Blitz & Spergel 1991, Stanek et al. 1994), the evidence from infrared surveys for a 
larger stellar bar in the inner disk (Weinberg 1992, Lopez-Corredoira et al. 1997, Benjamin 
et al. 2005), the discovery of the Sagittarius dwarf (Ibata et al. 1994; Ibata & Gilmore 
1995), and a significant asymmetry of unknown origin in the distribution of faint blue 
stars in Quadrant 1 (Ql) of the inner Galaxy (Larsen & Humphreys 1996). Each of these 
observations provides a significant clue to the history of the Milky Way. When combined 
with the growing evidence for Galactic mergers in addition to the Sagittarius dwarf, i.e. 
the Monoceros stream (Newberg et al. 2002, Ibata et al. 2003), the Canis Major merger 
remnant (Martin et al. 2004), the Virgo stream (Vivas et al. 2001, Martinez-Delgado et al. 
2005), we now realize that the structure and evolution of our Galaxy have been significantly 
altered by mergers with other systems. Indeed the population of the Galactic Halo and 
possibly the Thick Disk as well, may be dominated by mergers with smaller systems. 

The asymmetry in Ql, / = 20° — 45°, first recognized by Larsen & Humphreys (1996) 
was based on a comparison with complementary fields in the fourth Quadrant (Q4) using 
star counts from the Minnesota Automated Plate Scanner Catalog of the POSS I (Cabanela 
et al. 2003) 1 . To map the extent and shape of the asymmetric distribution and further 
identify the contributing stellar population, Parker et al. (2003) extended the search to 40 
contiguous fields from the digitized POSS I on each side of the Sun-Center line plus the 
same number of fields below the plane in Ql. They examined the star count ratios for 
paired fields in three color ranges: blue, intermediate and red. Over 6 million stars were 
used in the star count analysis covering almost 2000 square degrees on the sky. They found 



x The MAPS database is online at: http://aps.umn.edu 
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a 25% excess in the number of probable Thick Disk stars in Ql, I 20 to 60° and 20° to 
40° above and below the plane compared to the complementary fields (I 340 to 300°) 
in Q4. While the region of the asymmetry is somewhat irregular in shape, it is also fairly 
uniform and covers several hundred square degrees. It is therefore a major substructure in 
the Galaxy due to more than small scale dumpiness. Assuming that these are primarily 
main sequence thick disk stars, with a typical magnitude completeness limit of V ~ 18 mag, 
they are 1-2 kpc from the Sun and about 0.5 to 1.5 kpc above (and below) the plane. 

Parker et al. (2004) also found an associated kinematic signature. Using velocities 
from spectra for more than 700 stars, they not only found an asymmetric distribution in 
the vlsr velocities, but the Thick Disk stars in Ql have a much slower effective rotation 
rate uj compared to the corresponding Q4 stars. A solution for the radial and tangential 
components of the vlsr velocities reveals a significant lag of 80 to 90 km s _1 in the direction 
of Galactic rotation for the Thick Disk stars in Ql. 

Three possible explanations for the asymmetry are the fossil remnant of a merger, a 
triaxial Thick Disk or inner Halo, and interaction of the Thick Disk/inner Halo stars with 
the bar in the Disk. Given the lack of any spatial overlap with the path of the Sagittarius 
dwarf through the Halo (Ibata et al. 2003), and the predicted path of the Canis Major 
dwarf (Martin et al. 2004), its association with either feature is unlikely. Furthermore 
its spatial extent and apparent symmetry relative to the plane also do not automatically 
support a merger interpretation. Our line of sight to the asymmetry is interestingly in the 
same direction as the stellar bar in the Disk (Weinberg 1992; Benjamin et al. 2005), and the 
near end of the Galactic bar (Lopez-Corredoira et al. 1997; Hammersley et al. 2000), but 
the Disk bar is approximately 5 kpc from the Sun in this direction. Thus the stars showing 
the excess are between the Sun and the bar, not directly above it. However the maximum 
extent of the star count excess along our line of sight is not known. 
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Recently Larsen et al. (2008) showed the stellar density distribution from the MAPS 
scans for the Parker et al. (2003) fields in Ql and Q4 above the plane, demonstrating that 
the excess was in Ql and was not due to a ring above the plane (Juric et al. 2008). Larsen 
et al named the asymmetry feature the Hercules Thick Disk Cloud in recognition of the 
direction on the sky where the star count excess is maximum. Thus, interpretation of the 
Hercules Thick Disk Cloud is not clear-cut. While it might well be a fossil merger remnant, 
the star count excess is also consistent with a triaxial Thick Disk or inner Halo as well as a 
gravitational interaction with the stellar bar especially given the corresponding asymmetry 
in the kinematics. The distinction between a triaxial Thick Disk and interaction with a 
Disk bar may be difficult to discern. Indeed, it is unclear which may have formed first, and 
one or both may be the result of mergers. 

The release of the SDSS Data Release 5 (DR5) photometry in the direction of the 
observed asymmetry in Ql led to the discovery of a feature at much fainter magnitudes, the 
distant Hercules-Aquila cloud (Belokurov et al. 2007) and also confirmed the nearer star 
count excess in the inner Galaxy (Juric et al. 2008). The SDSS survey however is not well 
designed for a comprehensive survey of the Thick Disk inside the Solar orbit. It extends 
below b = 30° in only a few directions in Ql and has only limited coverage in Q4. It does 
not have the leverage in Galactic longitude needed to discriminate among the possible 
causes of the Hercules Thick Disk Cloud. 

To further explore the possible origin of the observed spatial and kinematic asymmetry, 
we have obtained multicolor CCD images to extend the star counts to fainter magnitudes 
to map the spatial extent of the asymmetry along the line of sight as a function of Galactic 
longitude and latitude. For example, if the Thick Disk is triaxial we would expect to 
observe the star count excess out to greater longitudes, but the star count excess described 
by Parker et al. (2003) appears to terminate near / ~ 55°. However, if it is triaxial in Ql 
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with its major axis above the bar, the stars would be further away at the higher longitudes. 
By extending the star counts to fainter magnitudes, corresponding to greater distances, we 
can search for the asymmetry at higher longitudes, from I of 50° to 75°. 

In this first paper we describe our observing program and present results on a search 
for triaxiality of the inner Halo and Thick Disk. Our second paper will discuss the star 
counts across the full range of longitudes in Ql and Q4 containing the asymmetry feature 
and a third paper in the series will analyze the kinematics of the associated stars. In the 
next section we describe the CCD observing program and data reduction. Determination of 
the star count ratios and their errors for the fields at the higher longitudes are discussed in 
§3. Our results do not support a triaxial interpretation of the asymmetry feature and are 
summarized in §4. 



2. Observations and Data Reduction 

To map the asymmetry feature both above and below the Galactic plane we identified 
67 fields ranging in longitude from I = 20° to 75°, and / = 340° to 285° and latitude 
b = ±20° to ±45°. The total sky coverage of the program is 47.5 square degrees. The 
distribution of our program fields on the sky is shown in Figure 1 and the information for 
each field is given in Table 1 with their Galactic and equatorial centers, instrument, field 
size and date of observation 2 . In this paper we will be discussing a subset of the fields 
relevant to the triaxiality question. These fields are circled on Figure 1 and are listed in 
Table 6. Our program required observations in both hemispheres. 

2 Each field is identified according to its Galactic longitude and latitude with an H to 
distinguish it from similar fields based on APS data 
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2.1. Northern Hemisphere - Steward Observatory Bok 2.3 Meter 

We used the Steward Observatory 2.3 meter Bok telescope with the 1 square degree 
90Prime camera in 2006 May, 2007 September and 2008 May. This mosaic camera has 4 
blue sensitive 4096x2096 CCDs. Counting cosmetic defects and the plate scale of 0.45" per 
pixel we had an effective field of view of approximately 1.02 square degrees. The CCDs are 
arranged in a square with intra-CCD gaps of about 10 arcminutes. We used a Johnson 
U,B,V + Cousins-Kron R filter set for the observations. This field of view was ideal for the 
project as it allowed each of the project fields to be observed in a single pointing. 

Program fields were observed twice in each filter with exposure times of 600 seconds 
in U, 480 seconds in B, and 90 seconds each in V and R. Exposure times were chosen to 
reach a limiting magnitude of 20th in U and 22nd in B, V and R. These exposures were of 
sufficient length to require autoguiding. The filters were not parfocal but focus was digital 
and corrections were easily applied. 

We reduced the data with the IRAF MSCRED mosaic image reduction package (Valdes 
1998). The data were corrected for cross talk, bias, and dark current. Flatfielding was done 
using night-sky flats from program fields supplemented by a series of random sky exposure 
with equivalent exposure times. CCD defects, bad pixels and satellite trails were removed. 

Astrometry was done with Larsen's MOSAF (Mosaic Astrometry Finder) program 
written for the Spacewatch 0.9 meter telescope (Larsen et al. 2007) which uses a 13 term 
polynomial for a plate model, the CFITSIO library (Pence 1999) and WCSLib from 
WCSTools (Mink 1999). The astrometric reference catalog was USNO-B1.0 (Monet et al. 
2003), and the residuals were approximately 0.3" in all frames. 

Photometric reductions used both Landolt Standards (Landolt 1992) and a secondary 
set of faint standards from Osmer et al. (1998), and followed the method of Hardie (1964). 
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The formal zero point errors were approximately 0.03 magnitudes per filter with the 
photometric RMS of 0.05 magnitudes per band. An example of the mean atmospheric 
extinction and color terms can be found in Table 2; seeing varied from 1.0" to 2.8" during 
the 2006 May observations. 

2.2. Southern Hemisphere - CTIO SMARTS 1.0 Meter 

Fields which could not be reached with the Bok 2.3 meter telescope were observed 
using the Y4KCam at the 1.0 meter SMARTS consortium telescope at CTIO in 2006 April 
and 2008 October. This camera is a 4096x4096 pixel single CCD with a plate scale of 
0.289" per pixel and a field of view of about 0.108 square degrees. Consequently 9 subfields 
were observed to yield an effective area of 0.95 square degrees with overlap. The Johnson 
U,B,V + Cousins-Kron R filter set was used for these observations. 

Program fields were observed once in each filter. During 2006 April, the exposure 
times were 300s in U, 150s in B, and 60s in the V and R bands. This was sufficient to reach 
19th magnitude in B. V and R and 18th magnitude in U. In 2008 October, the exposure 
times were increased to 400s in U, 485s in B, 160s in V and 190s in R to reach fainter 
magnitudes. These exposures were of sufficient length to require autoguiding. The filters 
were not parfocal but could be compensated by a focus reduction. 

Our initial processing of the frames was based on the IRAF scripts for the Y4KCam 
provided by Massey 3 . Flat fielding used both dome and twilight flats. Astrometry was 
done using WCSTools and the USNO-B1.0 catalog with results refined and checked against 
MOSAF. Seeing was moderate for the 2006 April observations, with FWHM from 1.5" to 
3.0" . 



3 Available at: http:/ /www. lowell.edu/users/massey/obins/y4kcamred. html 
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Photometric reductions used both Landolt Standards (Landolt 1992) and a secondary 
set of faint standards from Osmer et al. (1998). Photometric magnitudes were then 
determined using the same method as the 90Prime data to assure consistency. Examples of 
the mean atmospheric extinction and color terms are in Table 3. 

2.3. Catalog Creation 

Source detection, image parametrization and catalog creation used Source Extractor 
(Bertin & Arnouts 1996). After several tests to verify that we were detecting objects 
close to the image limits, the default size for the Gaussian convolution mask was used 
and magnitudes were determined using the MAG_AUTO feature which uses an adaptive 
aperture photometry routine. Aperture photometry was performed on the detected objects. 
Since most are stars, they were not subject to the same faint-end systematics as galaxies 
(Bernstein et al. 2002). We also discarded all objects detected within 4" of the edge of the 
frame as potentially truncated. 

A custom pipeline was written to apply the WCS astrometry solution and photometric 
calibration to each object in the catalogs for each field. We then used the WCS positions to 
match the objects observed in the different filters. To be included in the final catalog, each 
object had to be measured in at least two filters to provide color information. Problems 
with the star-galaxy classifier in Source Extractor are discussed in the next section. 

The final catalogs with magnitudes, positions, colors and object classifications will be 
available on-line with the second paper. 
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2.4. Star Galaxy Separation 

The neural network classifier in Source Extractor does not work on images with pixel 
scales similar to those in our data. We therefore developed our own simple classifier, called 
JAL, based on image parameters from Source Extractor sensitive to the differences between 
point sources and extended objects. We conducted a series of tests of the various image 
parameters to find three that were suitable for a simple parameter space cut which could 
be applied to each image. The following three parameter space comparisons were adopted: 

1. MU.MAX vs. IS O ARE A JM AGE captures the tendency of extended objects to 
have a surface brightness much lower than a star at the same isophotal diameter. A 
sample parameter space is shown in Figure 2, with the discriminator curve separating 
stars from galaxies. 

2. MAGJSO vs. MAG_APER isolates the tendency of extended objects to have a 
substantial fraction of its light outside of the core of the PSF, generating an isophotal 
magnitude larger than the aperture magnitude. An example of this parameter space 
is shown in Figure 3, with the discriminator curve separating stars from galaxies. 

3. XYJMAGE vs. MU.MAX shows how an extended object tends to have a 
moment in the intensity distribution much different than a star at the same peak 
surface brightness. The sample parameter space for these parameters is shown in 
Figure 4, with the discriminator curve separating stars from galaxies. This is the 
weakest of the three criteria and is therefore given one-half the weight of the other 
tests. 

We used these three discriminators in a simple Perl/Pgplot routine which allowed the 
user to interactively define a curve separating stars (scored as 0) from galaxies (scored as 
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1). The net classification for a given image was a weighted score of to 1 averaged from the 
three tests. During the final catalog matching an averaged net score was then computed 
from the V and R bands only. We conservatively considered an object with a average score 
greater than 0.3 to be a galaxy. 

We developed this method by comparison with other classifiers. Training and 
fine-tuning the parameter spaces was done using stars and galaxies classified with FOCAS 
in the Osmer and Hall standard fields observed each night. As a final test, we chose a 
90Prime field (H055+42) in the MAPS Catalog of the POSS I with object classification 
using a neural network (Odewahn et al. 1992, 1993) which was also observed and classified 
by SDSS. The results of our comparison as a function of magnitude are in Tables 4 and 
5 for the MAPS and SDSS, respectively. For V > 16, our classifier agreed with the APS 
classification 89% of the time or better and with SDSS 95% of the time or better. The 
lower agreement with the MAPS can be attributed to blended images on the MAPS scans 
which were classified as non-stellar or galaxies. The poor agreement with both APS and 
SDSS for V < 16 is easily explained by the onset of saturation in the 90Prime images for 
our exposure times. 

2.5. Comparison with SDSS 

Six of our program fields overlap SDSS. Figure 5 shows a photometry comparison 
between objects classified as stars by both SDSS and our classifier for the program field 
H055+42. The SDSS stars have been transformed to the standard V band (Rodgers et 
al. 2006). For this comparison we find a zero point difference of 0.02 magnitudes. For 
V < 20.0 the RMS difference between SDSS magnitudes and ours is 0.04 magnitudes, for 
20.0 < V < 21.5 the deviation is 0.09 magnitude and for V > 21.5 it is 0.31 magnitude. 
Our data is thus well-suited for our star count analysis down to V ~ 21. 



2.6. Extinction Corrections 



Our magnitudes and colors were corrected for interstellar extinction using the maps 
from Schlegel et al. (1998) and the standard interstellar extinction law (Cardelli et al. 
1989). Each star was corrected based on its coordinates using the interpolation routines 
provided. This approach assumes that the extinction comes from relatively near the Sun 
and is applied as a constant correction. Because our fields are intermediate to high galactic 
latitudes, extinction is relatively low. All but two fields have an average E(B — V) of less 
than 0.08. One high extinction field (H300-20) has an average E(B - V) of 0.36. The 
extinction-corrected color-magnitude diagram resembled the other fields with the placement 
of the blue ridge at the expected B — V color. Therefore extinction anomalies do not play 
a role in our evaluation of the star counts. 

2.7. Completeness Estimates 

The classical estimate of the completeness limit is determined from a plot of the log 
of the cumulative star counts vs. magnitude and determining the magnitude at which 
it deviates from the expected straight power law. For galactic structure work at the 
magnitudes we are considering however, this classical method does not work well because 
the density of Thick Disk stars changes with distance above the plane of the Galaxy. At 
fainter magnitudes there are simply fewer stars to be found at larger distances which can 
look like incompleteness. 

To derive our completeness estimates we therefore adopt a model-based estimate of 
completeness using Larsen's galactic model program GALMOD (Larsen & Humphreys 
2003). In Figure 6 we show the results of a completeness estimate for field H300-20 using 
a classical power law compared with one using the galactic model. While the power law 
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deviates by more than 10% from the observed cumulative counts as bright as V = 18.2, 
the model based estimate shows that due to the decreasing density of stars along the line 
of sight the completeness limit of the field is closer to V = 19.6. We use the model-based 
completeness estimate for the rest of our analysis. Our completeness limits are tabulated 
in Table 6. If the completeness limit appears fainter than V = 21.0, it is reported as being 
21.0. In general, due to their brighter limiting magnitudes, our CTIO fields are not as 
complete as the fields observed with the Bok telescope. 



3. The Star Counts Ratios - A Test for Triaxiality 

The fields used to test for triaxiality are circled on Figure 1 and listed in Table 6. For 
this initial analysis we use V and the B — V color only. A sample color-magnitude diagram 
for one of our fields is shown in Figure 7. 

To isolate the faint blue stars in the Thick Disk/inner Halo which show the stellar 
excess, we use the easily recognized "blue ridge" which identifies the locus of the main 
sequence stars in each field. The green lines in Figure 7 delineate this magnitude range 
over which all of the Y4KCam and 90Prime fields are complete. At these magnitudes, a 
well-defined peak in the number of stars is apparent near B — V ~ 0.6, the blue ridge, as 
illustrated in the color histogram in Figure 8. As we have discussed in previous papers 
(Larsen & Humphreys 1996; Parker et al. 2003), the GALMOD model shows that bluewards 
of this peak, we isolate a sample of inner halo and thick disk stars in the magnitude range 
16 < V < 19. Assuming that these are main sequence subdwarfs, these stars will be 
between 1 and 5 kpc from the Sun. 

We define the peak of the blue ridge, (B — V)p, as the median color of the stars within 
± 0.4 dex of the maximum color bin, see Figure 8. The resulting blue ridge peak colors are 
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included in Table 6. With the magnitude and color range thus defined, we also give the 
area of each field, the estimated completeness limit Vc, and the number of stars bluer than 
(B - V) P in the magnitude ranges 16 < V < 19, 17 < V < 18, and 18 < V < 19. 

We then determined the star count ratios for the paired fields in Ql and Q4 normalizing 
for the area of the field. The ratios for each magnitude range are summarized in Table 7 
with their error and the significance parameter, s computed as per Parker et al. (2003). 
The predicted ratios from GALMOD are also included. GALMOD ratios were required for 
all comparisons which involved fields above and below the plane because the Sun is not in 
the mid-plane of the disk but instead lies slightly above it (Humphreys & Larsen 1995). 
Ratios of fields above the plane compared to fields below the plane should be slightly less 
than one due to this geometry and the Sun's location in the Northern Galactic hemisphere. 
Comparisons between lines of sight in the same Galactic hemisphere should have a ratio 
near unity if the stellar distribution in the Galaxy is symmetric. 

There are four comparisons in Table 7: Ql and Q4 fields above the plane, Ql and Q4 
fields below the plane, Ql fields above with those below the plane, and the Q4 fields above 
and below the plane. 

We find a statistically significant excess in the ratios for the Ql to Q4 fields above the 
plane for the two fields at the lowest longitudes, I of 45° and 50° that are an extension of 
the Hercules Thick Disk Cloud asymmetry. The fields at the greater longitudes, I of 55°, 
60°, 65° and 75° show no trend of systematic deviations from a ratio of 1.0. For Ql above 
and below the plane, there is a small excess above the plane in the two lowest longitude 
fields that does not extend to higher longitudes. We do not find evidence for an excess when 
comparing Ql to Q4 below the plane or when comparing the Q4 fields above and below 
the plane. A few deviations from the expected ratios exist in isolated directions/magnitude 
bins and are of interest, but they do not show a systematic trend. 
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4. Summary and Conclusions — A Search for a Triaxial Thick Disk 

Our star count ratios summarized in Table 7 show a small excess of faint blue stars at I 
of 45° and 50° but not at the greater longitudes including the faintest magnitude intervals. 
These results support Parker et al's (2003) conclusion that the asymmetry feature in Ql 
extends to about 55°. We thus find no evidence that the Hercules Thick Disk Cloud and the 
star count excess is due to a triaxial Thick Disk. The ratios at I of 45° and 50° return to 1.0 
at the faintest magnitude bin (18 < V < 19), suggesting that we may be seeing through the 
asymmetry feature; a possible clue to its extent along our line of sight in this direction. At 
V ~ 19 however, these faint blue stars will be about 2 kpc above the plane, and may not 
be in the Thick Disk. 

We also note that the above vs. below the plane ratios in Ql and Q4 and the Q1/Q4 
ratio below the plane are consistent with the predictions of the Galactic model. There is no 
star count excess. This is not in agreement with Parker et al. (2003) who concluded that 
the excess exists both above and below the plane. This may be evidence that the Hercules 
Thick Disk Cloud is actually a debris stream. However, considering the small number of 
field pairs and the fact that we are apparently at the edge of the Cloud at these longitudes, 
we consider this a tentative result pending the analysis of the full set of fields. 

Paper II will include an analysis of the full dataset, mapping the stellar excess in 
Galactic longitude and latitude and along the line of sight to determine the full extent of 
the Cloud. Our third paper will include the spectroscopy and analysis of the kinematics of 
the stars showing the excess. 
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Table 1. Field Information for the Thick Disk Asymmetry Project. 



Field Name I b Mean RA (J2000) Mean Dec (J2000) Instrument Area Run 01 



H020+20 


20.63° 


18.79° 


l 7 h 17 m 12 s 


-01°46'31" 


90Prime 


1.04 


2006 


H020+32 


20.53° 


30.79° 


16 ft 35 m 28 s 


+04°09'44" 


90Prime 


1.04 


2006 


H020+47 


20.30° 


45.70° 


15M2 m 32 s 


+11°17'18" 


90Prime 


0.77 


2008 


H020-47 


21.02° 


-48.25° 


21 ft 36 m 27 s 


-28°00'21" 


Y4KCam 


0.80 


2008 


H023+40 


23.34° 


38.69° 


lQh 12 m Q1 s 


+09°59'47" 


90Prime 


0.77 


2008 


H025+40 


25.38° 


38.83° 


lQh U m 2 r 


+ 11°28'33" 


90Prime 


1.04 


2006 


H027+37 


27.31° 


35.72° 


16 /l 28 m 36 s 


+ 11°28'52" 


90Prime 


0.77 


2008 


H027-37 


28.03° 


-38.18° 


20 /l 59 m 56 s 


-20°22'32" 


Y4KCam 


0.80 


2008 


H027+40 


27.25° 


38.73° 


16H7 m 24 s 


+ 12°44'26" 


90Prime 


0.77 


2008 


H030+20 


30.51° 


18.89° 


17 fe 34 m 43 s 


+06°29'21" 


90Prime 


1.04 


2006 


H030-20 


30.90° 


-21.17° 


19 ft 58 m 42 s 


-11°28'38" 


90Prime 


0.77 


2007 


H033+40 


32.68° 


38.76° 


16 /l 24 m 25 s 


+ 16°37'47" 


90Prime 


0.77 


2008 


H033-40 


34.21° 


-41.28° 


21 fc 19 m 20 s 


-17°02'39" 


Y4KCam 


0.48 


2008 


H035+32 


35.34° 


30.89° 


16 fe 58 m 07 s 


+ 15°34'43" 


90Prime 


1.04 


2006 


H035-32 


36.06° 


-33.06° 


20 /l 50 m 41 s 


-12°23T5" 


90Prime 


1.04 


2006 


H042-20 


43.34° 


-41.08° 


21 /l 30 m 28 s 


-10°43'38" 


90Prime 


0.77 


2007 


H042+40 


42.09° 


38.99° 


16 ft 34 m 08 s 


+23°35'27" 


90Prime 


1.04 


2006 


H042-40 


43.30° 


-41.07° 


21 /l 30 m 24 s 


-10°44'47" 


90Prime 


0.77 


2007 


H044+40 


43.99° 


38.92° 


16 ft 36 m 18 s 


+24°58'46" 


90Prime 


0.77 


2008 


H045+20 


45.43° 


19.02° 


17 ft 58 m 57 s 


+ 19°18'40" 


90Prime 


1.04 


2006 


H045-20 


45.99° 


-20.96° 


20 /l 24 m 20 s 


+01°04'50" 


90Prime 


1.04 


2006 
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Table 1 — Continued 



Field Name / b Mean RA (J2000) Mean Dec (J2000) Instrument Area Run C 



H046+45 


45.08° 


43.79° 


16 h 16 m 17 s 


+27°02'01" 


90Prime 


0.77 


200 


H048+45 


47.27° 


43.97° 


16H7"W 


+28°36'35" 


90Prime 


0.77 


200 


H050+31 


50.23° 


30.09° 


17 ft 20 m 37 s 


+27°19'10" 


90Prime 


1.04 


200 


H050-31 


51.19° 


-31.87° 


2lHl m 17 s 


-00°38'31" 


90Prime 


1.04 


200 


H053+42 


52.79° 


41.05° 


16 ft 34 m 02 s 


+32°02'05" 


90Prime 


0.77 


200 


H055+42 


54.87° 


41.20° 


16 ft 34 m 35 s 


+33°36'18" 


90Prime 


1.04 


200 


H055-42 


56.55° 


-42.86° 


21 /l 57 m 14 s 


-03°16'54" 


90Prime 


0.59 


20C 


H060+20 


60.32° 


19.27° 


18 ft 21 m 40 s 


+32°25'12" 


90Prime 


1.04 


200 


H060-20 


61.06° 


-20.69° 


20 /l 54 m 40 s 


+ 12°59'39" 


90Prime 


1.04 


200 


H065+31 


65.08° 


30.36° 




+39°46'56" 


90Prime 


1.04 


200 


H065-31 


66.33° 


-31.70° 


21 /l 41 m 43 s 


+09°44'33" 


90Prime 


0.78 


20C 


H075+20 


75.28° 


19.54° 


18M5 m 52 s 


+45°44'02" 


90Prime 


1.04 


200 


H075-20 


76.14° 


-20.50° 


21 ft 32 m 53 s 


+23°49'24" 


90Prime 


0.77 


20C 


H285-20 


285.37° 


-20.11° 


08Ml m 31 s 


-72°10'17" 


Y4KCam 


0.70 


200 


H285+20 


286.13° 


19.82° 


ll ft 23 m 00 s 


-39°45'21" 


Y4KCam 


0.90 


200 


H295-31 


294.94° 


-31.35° 


04 /l 48 m 47 s 


-82°01'58" 


Y4KCam 


0.90 


200 


H295+31 


296.39° 


30.60° 


12M9 m 37 s 


-31°22'40" 


Y4KCam 


0.90 


200 


H300-20 


300.34° 


-20.41° 


ll ft 28 m 40 s 


-82°19'09" 


Y4KCam 


0.85 


200 


H300+20 


301.12° 


19.50° 


12 ft 36 m 33 s 


-42°47'15" 


Y4KCam 


0.90 


200 


H305-42 


304.75° 


-42.67° 


00 fe 27 m 26 s 


-74°56'07" 


Y4KCam 


0.90 


200 


H305+42 


306.66° 


41.38° 


12 ft 58"W 


-20°52'02" 


90Prime 


1.04 


200 
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Table 1 — Continued 



Field Name / b Mean RA (J2000) Mean Dec (J2000) Instrument Area Run C 



H307+32 


308.71° 


41.23° 


13 h 04 m 44 s 


-20°53'57" 


90Prime 


0.77 


200 


H310+31 


311.27° 


30.30° 


13H9 m 48 s 


-31°29'45" 


Y4KCam 


0.90 


200 


H312+45 


314.28° 


44.14° 


13 fc 20 m 13 a 


-17°27'33" 


90Prime 


0.77 


200 


H314+45 


315.77° 


44.09° 


13 ft 24 m 38 s 


-17°18'20" 


90Prime 


0.77 


200 


H315-20 


315.40° 


-20.79° 


16 ft 51 m 22 s 


-76°49'47" 


Y4KCam 


0.90 


200 


H315+20 


316.13° 


19.22° 


13 ft 52 m 43 s 


-41°26'05" 


Y4KCam 


0.90 


200 


H316+40 


317.52° 


39.04° 


13 fe 34 m 48 s 


-21°56'24" 


90Prime 


0.77 


200 


H318-40 


318.02° 


-41.07° 


22 fe 12"W 


-71°50'00" 


Y4KCam 


0.90 


200 


H318+40 


319.46° 


39.17° 


13 ft 40 m 57 s 


-21°25'46" 


90Prime 


1.04 


200 


H325-32 


325.23° 


-32.96° 


20H0 m 13 s 


-70°11'34" 


Y4KCam 


0.90 


200 


H325+32 


326.19° 


31.05° 


14H5 m 31 s 


-27°16'07" 


Y4KCam 


0.92 


200 


H327+40 


328.88° 


38.89° 


14H0 m 41 s 


-19°12'19" 


90Prime 


0.77 


200 


H330-20 


330.46° 


-21.01° 


18M0 m 53 s 


-64°13'34" 


Y4KCam 


0.90 


200 


H330+20 


330.97° 


18.98° 


14 /l 59 m 22 s 


-36°05'31" 


Y4KCam 


0.90 


200 


H333-37 


28.09° 


-37.78° 


20 /l 58 m 26 s 


-20°11'59" 


Y4KCam 


0.21 


200 


H333+37 


334.20° 


35.82° 


14 fe 32 m 34 s 


-20°03'01" 


90Prime 


0.77 


200 


H333+40 


334.30° 


38.81° 


14 ft 26 m 45 s 


-17°23'50" 


90Prime 


0.77 


200 


H335-20 


335.09° 


-41.24° 


21 fe 09 m 48 s 


-60°32'57" 


Y4KCam 


0.20 


200 


H335-40 


335.07° 


-41.19° 


21 fe 09 m 26 s 


-60°34'26" 


Y4KCam 


0.53 


200 


H335+40 


336.19° 


38.97° 


U h 31 m 47 s 


-16°32'02" 


90Prime 


1.04 


200 


H337+40 


338.24° 


38.76° 


U h 37 m 5A s 


-15°52'59" 


90Prime 


0.77 


200 
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Table 1 — Continued 



Field Name / b Mean RA (J2000) Mean Dec (J2000) Instrument Area Run C 



H340-20 


340.52° 


-21.09° 


18 h 36 m 21 a 


-55°23'38" 


Y4KCam 


0.90 


200 


H340+20 


340.86° 


18.91° 


15 ft 36 m 16 s 


-30°46'05" 


Y4KCam 


0.88 


200 


H340+32 


341.02° 


30.87° 


15 fc 04 m 29 a 


-21°10'01" 


90Prime 


1.04 


200 


H340+47 


341.29° 


45.75° 


14 fc 30 m 10 a 


-08°47'15" 


90Prime 


0.77 


200 
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Table 2: Mean extinction and color terms from the 90Prime camera on the Bok 2.3 meter 

telescope for M ay of 2006. 

Band or Color Atmospheric Extinction Color Terms 
V 0.142 ±0.024 -0.032 ±0.010 

BV 0.089 ±0.011 1.125 ±0.018 

VR 0.049 ±0.008 0.930 ± 0.010 

UB 0.232 ±0.009 0.957 ±0.021 



Table 3: Mean extinction and color terms from the Y4KCam camera on the SMARTS 1.0 

meter telescope for April of 2006. 

Band or Color Atmospheric Extinction Color Terms 
V 0.163 ±0.017 0.096 ±0.013 

BV 0.123 ±0.007 0.847 ±0.033 

VR 0.052 ±0.019 0.912 ±0.025 

UB 0.245 ±0.030 0.877 ±0.057 
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Table 4: Star-Galaxy Separation Comparison between the JAL classifier and the Automated 
Plate Scanner for field H055+42. The table presents the number of classified stars (S) and 
galaxies (G) for each method along with the number of classifications in agreement. The 
last columns express the overall agreement and agreement on only the stellar classification 
between methods. 



B Mag 


APS 


APS 


JAL 


JAL 


Agree 


Disagree 


Agree 


Disagree 


Overall 


Stellar 


Range 


S 


G 


S 


G 


G 


G 


S 


S 


Agreement 


Agreement 


15 - 16 


4 





3 


1 





1 


3 





75% 


75% 


16 - 17 


74 


2 


73 


3 


1 


2 


72 


1 


96% 


97% 


17 - 18 


377 


33 


388 


22 


10 


12 


365 


23 


91% 


96% 


18 - 19 


635 


82 


653 


64 


45 


19 


616 


37 


92% 


97% 


19 - 20 


432 


191 


499 


124 


110 


14 


415 


81 


84% 


96% 


20 - 21 


234 


241 


345 


130 


110 


19 


213 


131 


68% 


91% 


21 - 22 


137 


194 


264 


67 


55 


12 


123 


139 


54% 


89% 
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Table 5: Star-Galaxy Separation Comparison between the JAL classifier and the Sloan Digital 
Sky Survey for field H055+42. The table presents the number of classified stars (S) and 
galaxies (G) for each method along with the number of classifications in agreement. The 
last columns express the overall agreement and agreement on only the stellar classification 
between methods. 



B Mag 


SDSS 


SDSS 


JAL 


JAL 


Agree 


Disagree 


Agree 


Disagree 


Overall 


Stellar 


Range 


S 


G 


S 


G 


G 


G 


S 


S 


Agreement 


Agreement 


15 - 16 





4 


3 


1 


1 








3 


25% 


0% 


16 - 17 


63 


13 


73 


3 


3 





63 


10 


86% 


100% 


17 - 18 


384 


26 


388 


22 


10 


12 


372 


16 


93% 


96% 


18 - 19 


660 


57 


653 


64 


41 


23 


637 


16 


94% 


96% 


19 - 20 


498 


125 


499 


124 


108 


16 


482 


17 


94% 


96% 


20 - 21 


344 


131 


345 


130 


116 


14 


330 


15 


93% 


95% 


21 - 22 


248 


83 


264 


67 


64 


3 


245 


19 


93% 


98% 
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Table 6: Extinction Corrected Star Counts. 



r ieiu. in auie 
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V c a i 




- V )p 


iu <^ v <^ ry 
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1 / \ l 7 v. lO 


18^- 1/ ^ in 
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TTn/ic; on 
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Do 


1 1 7n 


/lOQ 
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4Do 
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1 
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20, 
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Dl 


KA 7 
04 


104 


ofin 
zDU 


Trnfin on 
HUDU-zU 


1 


,02 


21. 


,0 


n 


DD 


zDsy 


n/i n 
y4y 


C^OI c 

OZ1 


wnRnj_on 
riuou+zu 


1 


,00 


21. 


.0 


n 


70 


osnQ 

ZoUo 


1 m q 
lUlo 


001 


nnfic i Q1 
HUDO+ol 


1 


.01 


20. 


,5 


U. 


7n 


nQ7 
yo ( 


Q/1 o 
o4Z 


o04 


HUDO-ol 


0. 


.72 


21. 


,0 
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DD 


7/1 1 
(41 
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ofin 
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TTn7c;j_on 
liU ( O+zU 
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.02 


21. 
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.93 
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fi7n 
D III 


7Q K 


TTOfic; on 


0. 


,73 
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1 7Qfi 
1 1 oo 
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Doz 


lTOn^_i_Q1 





.93 


19, 
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n 
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OQQ 
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QSQ 
OOO 


nzyo-oi 





.93 


19. 


,2 


U. 


fin 

Dy 


nQQ 
yoo 


Ofi K 

Zoo 


40z 


TTQnn on 
lioUU-zU 





.93 


18. 


,5 


n 
U, 


Oo 


zOoO 


nno 
yyz 


/in/i c 
4y4 


H300+20 





.93 


19. 


,2 


0. 


.62 


2522 


932 


529 c 


H305+42 


1 


.02 


21. 


.0 


0. 


,64 


525 


170 


242 


H310+31 





,91 


19. 


,2 


0, 


,58 


870 


289 


399 


H315-20 





,93 


19. 


,4 


0, 


,58 


3176 


1052 


1536 


H315+20 





,93 


19. 


,5 


0, 


,58 


2880 


980 


1418 



"Estimated Completeness Limit for the Field in the V Band. 
b Median Location for the blue ridge Line. 

c Due to incompleteness of H300-20, last bin is only to V=18.5. 



Table 7. Star Count Ratios for fields with \l\ > 45 — 75° from the Galactic Center. 



Field 


h/h 


61/62 


GALMOD Ratio Predictions 






Observed Count Ratios 






Ratio 




16 < V < 19 


17 < V < 18 


18 < V < 19 


16 < V < 19 


s 


17 < V < 18 


s 


18 < V < 19 


s 


Quadrant 1/Quadrant 4 ratios above the Galactic Plane 


H045+20/H315+20 


45/315 


+20/+20 


1.00 


1.00 


1.00 


1.25 ± U.Uo 


8.12 


1.6b ± U.Uo 


6.45 


1.U7 ± U.U4 


1.92 


H050+31/H310+31 


50/310 


+31/+31 


1.00 


1.00 


1.00 


1.21 ± U.Uo 


3.87 


1 q 1 _i_ n 1 n 
l.ol ± U.1U 


3.07 


1 _L n (17 

l.Uo ± U.U( 


0.65 


H055+42/H305+42 


55/305 


+42/+42 


1.00 


1.00 


1.00 


1.U4 ± U.Ub 


0.66 


n n 1 _i_ n 1 n 
u.yi ± U.1U 


0.90 


1 (17 _L n 1 n 
l.U ( ± U.1U 


0.70 


H060+20/H300+20 


60/300 


+20/+20 


1.00 


1.00 


1.00 


1 nn i n no 

l.Oo ± O.Oo 


1.19 


1.01 ± 0.05 


0.24 


0.97 ± O.Oo 


0.53 


H065+31/H295+31 


65/295 


+31/+31 


1.00 


1.00 


1.00 


0.97 ± 0.05 


0.67 


1.05 ± 0.08 


0.64 


0.85 ± 0.06 


2.37 


H075+20/H285+20 


75/285 


+20/+20 


1.00 


1 nn 
l.UU 


l.UU 


1.U1 ± U.Uo 


U.oU 


n c\i J_ n nc; 
u.y / ± U.Uo 


U.o / 


n nft J_ n nc; 
u.yo ± U.UO 


n Q1 
U.J51 


Quadrant 1/Quadrant 4 ratios below the Galactic 


Plane 


















H045-20/H315-20 


45/315 


-20/-20 


1.00 


l.UU 


i on 
l.UU 


i nn -U n no 

l.UU it U.Uz 


U.UD 


1 nc 4- n nc; 
l.Uo it U.Uo 


l.Do 


n on 4- n m 
u.yu it U.Uo 


o.lo 


H060-20/H300-20 


60/300 


-20/-20 


1.00 


1.00 


1.00 


0.97 ±0.03 


1.23 


0.87 ±0.04 


3.23 


0.96 ± 0.06 


0.64 


H065-31/H295-31 


65/295 


-31/-31 


1.00 


1.00 


1.00 


1.03 ± 0.03 


0.89 


1.02 ±0.06 


0.37 


1.02 ± 0.06 


0.43 


H075-20/H285-20 


75/285 


-20/-20 


1.00 


1.00 


1.00 


1.03 ± 0.05 


0.51 


1.17 ±0.10 


1.73 


0.77 ±0.06 


3.91 


Quadrant 1 ratios above/below the Galactic Plane 


















H045+20/H045-20 


45/45 


+20/-20 


0.94 


0.94 


0.97 


1.13 ± 0.03 


6.33 


1.18 ±0.05 


3.96 


1.10 ± 0.04 


3.25 


H050+31/H050-31 


50/50 


+31/-31 


0.97 


0.97 


0.98 


1.09 ± 0.05 


2.40 


1.15 ±0.08 


2.25 


1.01 ± 0.07 


0.43 


H060+20/H060-20 


60/60 


+20/-20 


0.94 


0.93 


0.96 


0.99 ± 0.03 


1.66 


1.06 ±0.05 


1.80 


1.07 ±0.07 


1.85 


H075+20/H075-20 


75/75 


+20/-20 


0.94 


0.93 


0.96 


0.84 ±0.05 


2.00 


0.84 ±0.05 


1.80 


0.86 ± 0.04 


2.50 


H065+31/H065-31 


65/65 


+31/-31 


0.97 


0.96 


0.98 


0.90 ±0.04 


1.75 


0.94 ±0.08 


0.25 


0.94 ±0.08 


0.50 


Quadrant 4 ratios above/below the Galactic Plane 


















H315+20/H315-20 


315/315 


+20/-20 


0.94 


0.94 


0.97 


0.91 ± 0.02 


1.50 


0.93 ±0.04 


0.25 


0.92 ± 0.03 


1.60 


H300+20/H300-20 


300/300 


+20/-20 


0.94 


0.93 


0.96 


0.99 ± 0.03 


1.60 


0.94 ±0.04 


0.25 


1.07 ±0.07 


1.56 


H285+20/H285-20 


285/285 


+20/-20 


0.94 


0.93 


0.96 


0.85 ± 0.04 


2.25 


0.88 ±0.05 


1.00 


0.91 ± 0.05 


1.00 



Table 7— Continued 



Field 
Ratio 


h/h 


bi/b 2 


GALMOD Ratio Predictions 
16 < V < 19 17 < V < 18 18 < V < 19 


16 < V < 19 


Observed Count Ratios 
s 17<V<18 s 


18 < V < 19 s 


H295+31/H295-31 


295/295 


+31/-31 


0.97 0.96 0.98 


0.95 ±0.04 


1.03 1.05 ±0.09 1.00 


0.85 ±0.06 2.16 
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Distribution of Triaxia I ity Fields (Aitoff Projection) 




Galactic Longitude I (degrees) 

Fig. I. — Distribution on the sky of project fields. Fields observed with 90Prime are plotted 
as open squared and fields observed with Y4KCam are represented as filled squares. Fields 
used for our first result on triaxiality are circled. The size of each icon is roughly twice as 
large as the actual sky coverage of the corresponding field. 
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First Star Galaxy Discrimination Criteria 




2000 4000 0000 8000 

SExtractor ISOAREAJMAGE 

Fig. 2.— Star/galaxy separation with MILMAX vs. ISOIMAGE_AREA. The stellar locus 
is well defined. The plot represents all points on amplifier "A" of 90Prime during a 120 
second R band exposure on the night of May 26, 2006. The discrimination line is plotted in 
green. Objects which are classified as a galaxy when all three discriminators are applied are 
plotted in red while objects which remain on the stellar locus are plotted in blue. 
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Second Star Galaxy Discrimination Criteria 
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Fig. 3. — Star/galaxy separation with MAG_ISO versus MAG_APER. The plot represents 
all points on amplifier "A" of 90Prime during a 120 second R band exposure on the night 
of May 26, 2006. The discrimination line is plotted in green. Objects which are classified 
as a galaxy when all three discriminators are applied are plotted in red while objects which 
remain on the stellar locus are plotted in blue. 
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Third Star Galaxy Discrimination Criteria 
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SExtractor MILMAX 

Fig. 4. — Star/galaxy with XYJMAGE versus MLLMAX. The plot represents all points on 
amplifier "A" of 90Prime during a 120 second R band exposure on the night of May 26, 
2006. The discrimination line is plotted in green. Objects which are classified as a galaxy 
when all three discriminators are applied are plotted in red while objects which remain on 
the stellar locus are plotted in blue. The spike in the XYJMAGE direction which occurs at 
MILMAX ~ 13.7 is due to charge bleed for the brightest objects. Since these objects are 
already saturated they are not contaminating our star count samples. 
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90 Prime V Magnitude 



Fig. 5. — Photometric comparison for H0055+42 between the 90Prime V band data and the 
SDSS DR5, converted to V band using the transformations of Rodgers et al. (2006). 
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Completeness Estimate, Integrated Counts for H300 — 20 
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Fig. 6. — Integrated counts for field H300-20 as observed from CTIO (blue curve). Magnitude 
limits of interest in this paper (16 < V < 19) are delineated with green lines. Completeness 
estimates are also plotted: a classic power law matched to the counts for V < 19 (magenta 
curve) and a model-based estimate using the best fit parameters from Larsen & Humphreys 
(2003) (orange curve). The power law implies a completeness limit more than a magnitude 
brighter the model-based method. 
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B-V.V Color Magnitude Diagram, H045 + 20 
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Fig. 7. — Example Color-Magnitude diagram in V and B-V for H045+20, observed using 
90Prime. The magnitude range used in this analysis is bounded in green and the location of 
the "blue ridge" is shown in blue. 
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B-V Color Histogram, H045 + 20 
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Fig. 8. — The number of stars vs. color for H045+20 for all stars with 16 < V < 19 The 
location of the blue ridge line as determined by our median routine is indicated. 



